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The present status of calculations of the nuclear matrix elements for neu- 
trinoless double beta decay is reviewed. A proposal which allows in principle 
to measure the neutrinoless double beta decay Fermi matrix element is briefly 
described. 



1. Introduction 

Neutrino oscillation experiments have proven that neutrinos are massive particles 
(see, e.g., Ref. PQ). However, the absolute scale of the neutrino masses cannot in 
principle be deduced from the observed oscillations. To determine the absolute 
neutrino masses down to the level of tens of meV, study of the neutrinoless double 
beta {Qvf3j3) decay ^El^ — > z+2 El/v-2 + 2e~ becomes indispensable. Furthermore, 
this process, which violates the total lepton number by two units, is an experimentum 
crucis to reveal the Majorana nature of neutrinos, i.e. if neutrino is identical with 
antineutrino (see, e.g., Ref. [21 EH HI [2] ) • 

Determination of the effective Majorana mass (or relevant GUT and SUSY parame- 
ters depending on what mechanism of the Ovflfl decay dominates) from experimental 
data on the 0z//3/3-decay lifetimes can be only as good as the knowledge of the nu- 
clear matrix elements M 0u on which the Qvf3f3 decay rates depend. Thus, a better 
understanding of the nuclear structure effects important for describing the matrix 
elements is needed to interpret the future data accurately. It is crucial in this con- 
nection to develop theoretical methods capable of reliably evaluating the nuclear 
matrix elements, and to realistically assess their uncertainties. 

In general, barring contributions different from light Majorana neutrino exchange, 
the inverse lifetime in a given nucleus is the product of three factors, 

{ify)- 1 = CP \M<»\ 3 m*,, (1) 

where G 0u is a calculable phase space factor, M 0l/ is the 0u/3{3 nuclear matrix element, 
and m/3/3 is the so-called "effective Majorana neutrino mass" which, in standard 
notation [6], reads mpp = |Xa=i m « 5 with m 8 and U e i being the neutrino masses 
and the v e mixing matrix elements, respectively. 

Two-neutrino double beta (2i , (3f3) decay, ^El^r — ► Z _^E\ N _ 2 + 2e~ + 2V e is a 
second-order weak process which probes the same mother and daughter nuclei as 



Ovf3f3 decay. It is allowed within the standard model and has been observed in sev- 
eral nuclei (see, e.g., [7]). The decay provides a particularly important benchmark. 
As it was extensively demonstrated in [8], the spread of QRPA calculation results 
for M Qu can be significantly reduced by constraining the nuclear model with the 
corresponding experimental 2v(3/3 decay lifetimes. 

The calculation of the matrix element M Qv for a candidate Qv/3/3 nucleus is notori- 
ously difficult. It requires the detailed description of a second-order weak decay from 
a double-even mother nucleus (Z, A) to a double-even daughter nucleus [Z + 2, A) 
via virtual states (with any multipolarity J n ) of the intermediate nucleus (Z + l, A). 
Two basic methods are used in the evaluation of M 0u , the quasiparticle random 
phase approximation (QRPA), with its various modifications 0, [9], [TUl [HI [12] and 
the nuclear shell model (NSM) [13] (very recently M 0v have also been calculated 
within the IBM [2]). The NSM aims at complete describing the nuclear wave 
functions by taking into account nucleon configurations of all possible complexity. 
By diagonalizing the nuclear Hamiltonian within the model space the energies of 
ground and excited states of nucleus as well as the corresponding wave functions can 
be calculated. However, application of the NSM to description of the medium and 
heavy nuclei beyond the p/-shell immediately faces factorial growth of the dimension 
of the model space . Therefore, severe truncation of the size of the single-particle 
basis is usually made in the medium and heavy nuclei. 

In contrast to the NSM, within the QRPA and the renormalized QRPA (RQRPA) 
one can include essentially unlimited set of single-particle states, but only a limited 
subset of configurations (iterations of the particle-hole, respectively two-quasiparticle 
configurations) is taken into account. On the other hand, within the QRPA there 
is no obvious procedure that determines how many single particle states one should 
include. Hence, various authors choose this crucial number ad hoc, basically for 
reasons of convenience. 



2. Current status of calculations of M Uz/ 

The Ov/3/3 nuclear matrix element is a sum of Fermi (F), Gamow- Teller (GT) transi- 
tions and (a small) tensor (T) contributions (see, e.g., [8J for detailed representations 
of the contributions), 

M 0u = M°" T + M° u 1- , (2) 

9a 

and nuclear models are needed to estimate the different components (X = F, 
GT, T). In the above expression, g& is the effective axial coupling in nuclear matter, 
not necessarily equal to its "bare" free-nucleon value g& — 1.25. A direct comparison 
of M 0v calculated for different qa can be done in terms of the matrix element M' 0u = 

There has been great progress in the calculations of M 0u over the last five years. 
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Figure 1: Current status of calculations of M 0u for the light neutrino exchange 
mechanism (with the Jastrow-like s.r.c. and gA = 1.25) within different nuclear 
structure models (QRPA: [8] (T'07) and [JJ (J'07); NSM [TJ]; IBM2 [TJ]). 



A comparison of the results by different groups (with the Jastrow-like short-range 
correlations (s.r.c.) and with the unquenched value = 1.25) is represented in 
Fig. [TJ One can see in the figure that the matrix elements M 0v of different groups 
calculated within the QRPA seem to converge. At the same time, the M 0u of the 
NSM are systematically and substantially smaller (up to a factor of 2 for lighter 
nuclei like 76 Ge) than the corresponding QRPA ones. There is now an active dis- 
cussion in literature on what could be the reason of such a discrepancy, a too small 
single-particle model space of the NSM or a neglect of complex nuclear configura- 
tions within the QRPA. The recent results obtained within the IBM method [H] 
agree surprisingly well with the QRPA ones of Ref. [8]. 

As already mentioned, all models employ truncations: the NSM severely truncates 
the s.p. model space whereas the QRPA does so with respect to the configuration 
space. The question relevant for calculation of the Qv/3/3 nuclear matrix elements is 
which truncation induces the smallest error in M 0u . 

In Ref. [13] the difference between the NSM and the QRPA results is attributed to a 
neglect of a subset of ground state correlation in the latter. The claim as it appears 
in [TJ] is "the QRPA can be said to be a "low seniority approximation", roughly 
equivalent to the s < 4 ISM truncations, that overestimate the NME's However, 



it is obviously incorrect that all the components with the seniority s > A are not 
included within the QRPA since an analytic expression for the QRPA ground state 
as a coherent state built on top of the BCS vacuum can explicitly be derived (see, 
e.g., [E]) This representation shows that all the configurations with the seniority 
AN (where N = 0,1,2,...) are taken into consideration when the QRPA ground 
state is calculated. Of course, those represent only a part of the entire ground state 
correlations but the part which is the most relevant for calculating the transition 
amplitudes. 

From our point of view it is namely too small single-particle model space used in 
the NSM calculations that is responsible for the suppression of the calculated NSM 
M 0u . Only very limited number of negative parity configurations can be constructed 
within the Ohu model space. Therefore, contribution from many 0^/3/3-transitions 
(dipole, spin-dipole etc.) via negative parity intermediate states are missing in the 
NSM description. Such transitions contribute a lot to the M 0u as demonstrated by 
the QRPA results. Thus, it is natural to expect the NSM matrix elements M 0l/ 
to come out small because some important transitions contributing to the matrix 
element simply cannot be described in such a small basis. 

3. QRPA analysis of uncertainties in 

At present, the most elaborate analysis of uncertainties in the Ov/3/3 decay nuclear 
matrix elements M Qu calculated within the QRPA and the RQRPA has been per- 
formed in Refs. [9|, [10]. Single-particle model spaces comprising iV =2, 3 and 5 
major oscillator shells were used in the (R)QRPA calculations along with different 
representations of the short range correlations. The experimental 2uf5f3 decay rates 
were used to adjust the most relevant parameter, the strength g pp of the particle- 
particle interaction, and thus to "calibrate" the QRPA estimates of M 0u . The major 
observation of Refs. [HJ [9j [10] is that such a procedure makes the calculated M 0u 
essentially independent of the size of the single-particle basis of the QRPA. Further- 
more, the matrix elements have been demonstrated to also become rather stable 
with respect to the possible quenching of the axial vector coupling constant g^. 
Despite the fact that the 2vf}(3 decay process probes only a subset of the intermediate 
states relevant for OvfiP decay (i.e., only those with J n = 1 + , via GT transitions), 
it is just the 1 + contribution to the total Ovflfi matrix element that reveals a pro- 
nounced sensitivity to g pp , in contrast to the other multipole contributions [8]. This 
observation justifies the aforementioned fitting procedure employed in Ref. [8]. 
The matrix elements M' 0u calculated for the three single-particle bases and a fixed 
gA are relatively close to each other. For each nucleus the corresponding average 
(M /0iy ) matrix elements (averaged over the three choices of the single-particle space) 
was evaluated in Ref. [8], as well as its variance a. The final (R)QRPA results 
obtained by using unquenched and quenched values of g^ (gA = 1-25 and g^ = 1-0, 
respectively) are presented in graphical form in Fig. [2J The full uncertainties of the 
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Figure 2: Average nuclear matrix elements (M' 0u ) and their variance (including the 
error coming from the experimental uncertainty in M 2u ) for both the QRPA and 
the RQRPA M. 



calculated M' v includes also the ones induced by the experimental uncertainties in 
M 2v . One can see that not only is the variance substantially less than the average 
value, but the results of QRPA are quite close to the RQRPA values. 
A systematic analysis of the effect of different choices of the s.r.c. on M 0u calcu- 
lated within the QRPA and the RQRPA has been performed in Refs. [SJ [9], [TO] . 
The QRPA results obtained with the Jastrow-like and the unitary correlation op- 
erator method (UCOM) treatment of the s.r.c. are represented in Fig. [31 Also 
shown in the figure are the results of the first self-consistent calculation [TUj which 
uses the residual nuclear interactions as well as the s.r.c. derived from the same 
modern realistic nucleon-nucleon potentials, namely from charge- dependent Bonn 
potential (CD-Bonn) and the Argonne V18 potential. Larger matrix elements M 0u 
are obtained as compared with the traditional approach of using the Jastrow-like 
treatment of the s.r.c. 
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Figure 3: Effect of different choices of the s.r.c. on M 0l/ calculated within the 
QRPA El HD! (g A = 1-25). 

4. Can one measure M 0zy ? 

Given the uncertainty in calculated M 0u of different groups, it is of great importance 
if there are experimental means allowing to measure these matrix elements. In a 
recent work Ref. [TH] a proposal is put forward which shows that such a measurement 
of the Fermi part Mp" of the total matrix element M 0v is in principle possible. Here, 
we would like to briefly discuss this result. 

A similarity between the 0vf3f3 neutrino potential and the radial dependence of the 
two-body Coulomb interaction is exploited in Ref. [16] . When in addition one makes 
use of the isospin conservation by strong interaction, the matrix element Mp u can be 
transformed as to acquire the form of an energy-weighted double Fermi transition 
matrix element which is dominated by the amplitude of the double Fermi transition 
via the IAS in the intermediate nucleus into the ground state of the final nucleus: 



--u IAS {0 f \f-\IAS){IAS\f-\0 i ). 



(3) 



Here, the second Fermi transition amplitude is due to an admixture of the double IAS 
in the final nucleus to the ground state of the parent nucleus: (Of\T~\IAS) (IAS\T~\0 i ) 
(O f \DIAS)(DIAS\f-\IAS)(IAS\f-\Oi). 



Therefore, the total M% v can be reconstructed according to Eq. ([3]), if one is able 
to measure the AT = 2 isospin-forbidden m.e. (Of\T~\IAS), for instance in charge- 
exchange reactions of the (n,p)-type. Using recent QRPA calculation results for 
Mp , this m.e. can roughly be estimated as (Of\T~\IAS) ~ 0.005, i.e. about a 
thousand times smaller than the first-leg m.e. (IAS\T~\0i) « — Z. This strong 
suppression of (0/|T~| IAS) reflects smallness of the isospin mixing effects in nuclei. 
The IAS has been observed as a prominent and extremely narrow resonance and its 
various features have well been studied by means of (p,n), ( 3 He,t) and other charge- 
exchange reactions. This gives us a hope that a measurement of (0f\T~\IAS) in the 
(n,p) charge-exchange channel might be possible. More generally, a measurement 
by whichever experimental mean of the AT = 2 admixture of the DIAS in the final 
ground state would be enough to determine M F V . 

Of course, by measuring only M F U one does not get the total m.e. M 0u but rather 
its sub-leading contribution. However, knowledge of M F V itself brings a very impor- 
tant piece of information, since it can help to discriminate between different nuclear 
structure models in which calculated M F V may differ by as much as a factor of 5. 
In addition, the ratio Mf '/M^ T may be more reliably calculable in different models 
than M F U and Mq^ separately. Let us put forward here some simple arguments 
in support of the latter statement. Since only small internucleon distances deter- 
mine M 0u , then only nucleon pairs in the spatial relative s-wave must dominantly 
contribute to the m.e.. The isotensor Coulomb interaction only couples T = 1 
pairs which must then be in the state with the total spin S = to assure anti- 
symmetry of the total two-body wave function. Because of this and the fact that 
<Ti • cr 2 \S = 0, T = 1) = —SIS' = 0, T = 1), a natural estimate for the Gamow- Teller 
m.e. is Mq? t = — 3M F U provided the neutrino potential is the same in both F and 
GT cases. The high-order terms of the nucleon weak current which are present in 
the case of the GT m.e., but absent in the F m.e., change a bit this simple estimate 
to Mq , t /M f v ~ —2.5. Also, an uncertainty of few per cents may come from the 
difference in the mean nuclear excitation energies in the F and GT cases. It is worth 
noting that the recent QRPA results [H [HI [12] are in good correspondence with 
these simple estimates. 



5. Conclusion 

In this contribution the present status of calculations of the nuclear matrix elements 
M 0v for neutrinoless double beta decay is reviewed. The matrix elements M 0u of 
different groups calculated within the QRPA seem to converge. At the same time, 
the M Qv of the SM are substantially smaller (up to factor 2 for lighter nuclei like 
76 Ge) and this discrepancy is under active discussion in literature now. The recent 
results obtained within the IBM method agree surprisingly well with the QRPA 
ones. In addition, we have briefly described a proposal which allows in principle to 



measure the neutrinoless double beta decay Fermi matrix element. 
Acknowledgement: Many of the original results reported here have been obtained 
in a long-time fruitful collaboration with Amand Faessler, Fedor Simkovic and Petr 
Vogel, which is gratefully acknowledged by the author. The work is supported in 
part by the DFG within the SFB TR27 "Neutrinos and Beyond" . 



[1] R. D. McKeown and P. Vogel, Phys. Rep. 394, 315 (2004). 

[2] F. Boehm and P. Vogel, Physics of Massive Neutrinos 2nd ed., Cambridge Univ. 
Press, Cambridge, 1992. 

[3] A. Faessler and F. Simkovic, J. Phys. G 24, 2139 (1998). 

[4] J. Suhonen and O. Civitarese, Phys. Rep. 300, 123 (1998). 

[5] Frank T. Avignone III, Steven R. Elliott, and Jonathan Engel, Rev. Mod. Phys. 
80, 481 (2008). 

[6] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 
[7] A.S. Barabash, Czech. J. Phys. 50, 437 (2006). 

[8] V.A. Rodin, A. Faessler, F. Simkovic and P. Vogel, Phys. Rev. C 68, 044302 
(2003); V. A. Rodin, A. Faessler, F. Simkovic and P. Vogel, Nucl. Phys. A766, 
107 (2006); ibid. A793, 213(E) (2007); 

[9] F. Simkovic, A. Faessler, V.A. Rodin, P. Vogel, and J. Engel, Phys. Rev. C 77, 
045503 (2008); 

[10] F. Simkovic, A. Faessler, H. Muther, V.A. Rodin, and M. Stauf, Phys. Rev. C 
79, 055501 (2009). 

[11] V. Rodin and A. Faessler, Phys. Rev. C 77, 025502 (2008). 

[12] M. Kortelainen, O. Civitarese, J. Suhonen, and J. Toivanen, Phys. Lett. B647, 
128 (2007); M. Kortelainen and J. Suhonen, Phys. Rev. C 75, 051303(R) (2007); 
M. Kortelainen and J. Suhonen, Phys. Rev. C 76, 024315 (2007). 

[13] E. Caurier, F. Nowacki, A. Poves, Eur. Phys. J. A 36, 195 (2008); E. Caurier, 
J. Menendez, F. Nowacki, A. Poves, Phys. Rev. Lett. 100 052503 (2008); J. 
Menendez, A. Poves, E. Caurier, F. Nowacki, Nucl. Phys. A 818, 139 (2009). 

[14] J. Barea and F. IacheUo, Phys. Rev. C 79, 044301 (2009). 

[15] P. Ring and P. Schuck, The Nuclear Many Body Problem ( Springer- Verlag, 
Berlin, 1980). 

[16] Vadim Rodin and Amand Faessler, arXiv:0906.1759 [nucl-th]. 



